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The structural description of disordered systems has been a longstanding challenge in physical science. We 
propose an atomic cluster alignment method to reveal the development of three-dimensional topological 
ordering in a metallic liquid as it undercools to form a glass. By analyzing molecular dynamic (MD) 
simulation trajectories of a Cu 6 4. 5 Zr 3 5 5 alloy, we show that medium-range order (MRO) develops in the 
liquid as it approaches the glass transition. Specifically, around Cu sites, we observe "Bergman 
triacontahedron" packing (icosahedron, dodecahedron and icosahedron) that extends out to the fourth 
shell, forming an interpenetrating backbone network in the glass. The discovery of Bergman-type MRO 
from our order-mining technique provides unique insights into the topological ordering near the glass 
transition and the relationship between metallic glasses and quasicrystals. 

A correct description of the order parameter in a liquid or glass system is the key to understanding and 
predicting the dynamics of glass formation and phase selection in undercooled liquid materials. However, 
the lack of long range order (LRO) in the system introduces a spatially varying rotation between the 
structural motifs found at different positions in the system making it difficult to discern the details of the 
developing structural order both in scattering experiments as well as in atomistic simulations 1 " 14 . The directional 
averaging associated with the measurement of simple radial correlations in scattering experiment obscures the 
angular and species-specific information necessary to quantify the topological and chemical order even in the first 
coordination shell 1 " 5 . 

In the analysis of atomistic liquid and glass geometries, although several widely used methods (e.g. Bond- 
Orientation order 15 , Honeycutt- Anderson index 1617 and Voronoi tessellation 18 ) exist to complement the struc- 
tural information provided by radial-type pair correlations, they suffer from several common deficiencies: (i) 
They do not provide a visual picture of the structural motifs of the system although some features can be inferred 
from the numerical statistics provided by these analyses; (ii) They only address structural arrangements in the first 
coordination shell, extension to characterize medium range order beyond the first shell is not straightforward or, 
in many cases, non-existent; (iii) They do not treat different structural motifs on an equal footing, which can make 
comparison between different structural motifs difficult. For example, Voronoi indices can describe icosahedral 
short range order (ISRO) very well but run into mathematical difficulties when describing fcc-type short range 
order (SRO); (iv) They provide population statistics but do not provide a quantitative measure of the distance of a 
given cluster in the system from different competing SRO templates. 

In this paper, we propose a robust cluster alignment tool that corrects all of the above deficiencies in existing 
methods. Our method provides a visual picture of the structural motif as well as a scheme for defining an order 
parameter to measure the strength of any structural motif at a given locality in the system. The definition of this 
quantitative order parameter treats all different structural motifs on an equal footing and can be incorporated into 
modified MD schemes to manipulate the tendencies of the system towards various competing structural motifs. 
By applying this alignment method to the well-studied Cu 6 4. 5 Zr 3 5 5 alloy where the MRO still remains unclear 819 " 21 , 
we clearly demonstrate the development of Bergman triacontahedron MRO upon cooling from the liquid state. 
The discovery of Bergman-type MRO from our order-mining technique provides unique insights into the nature 
of the glass transition and the relationship between metallic glasses and quasicrystals. Our "order mining" method 
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Figure 1 | Atomic distribution of the Cu-centered clusters from the MD simulation of Cu 64 . 5 Zr35.5 at 300K. (a) as-simulated atom positions for 2,000 
randomly selected Cu-centered clusters each contains 301 atoms; (b) atom positions after the alignment. For clarity the image is a 2 A thick slice through 
the center of the clusters with Cu atoms in blue and Zr atoms in yellow. 



is generic and can be used to uncover hidden order in any form of 
disordered systems including metallic liquids and glasses. 

Results 

The 3D atomistic structures for both the Cu 64 _ 5 Zy 35 5 liquid and glass 
are generated by MD simulations using an accurate interatomic 
potential based on the Finnis- Sinclair method (FS) 21 . Each structure 
contains 5000 atoms in a cubic simulation cell, with periodic bound- 
ary conditions. The glass structures are generated by cooling from the 
equilibrated high temperature liquid at essentially zero pressure 
using a cooling rate as low as 4.91 X10 9 K/s. Details about the 
MD simulations and comparison with scattering experiments are 
reported elsewhere 22 . Structures at five temperatures (1500K, 
HOOK, 900K, 700K and 300K) are used in the present analysis, where 
the structures at 700K and 300K are below the glass transition tem- 
perature which is around 750K 21 . 

In order to reveal the development of three-dimensional topo- 
logical ordering upon cooling in Cu 64 5 Zr 35 5 liquid and glass, 2000 
atomic clusters each of a Cu atom in the center and about 300 atoms 
around the center atom are randomly picked up from the MD simu- 
lation trajectory at each temperature. The 2000 Cu-centered clusters 
at a given temperature are then aligned with each other to minimize 
the overall mean square distances among the atoms in different 
clusters so that the clusters will rotate to have a common orientation. 
The resulting alignment reveals the inherent structural motif char- 
acterizing the local environment of the given chemical constituent. 
More details about the alignment procedure are given in the methods 
section and Ref. 22 . 



An example of atomic distribution before and after the alignment 
is shown in Fig. 1. In this plot, 2000 Cu centered atomic clusters, each 
consisting of their 300 neighboring atoms, are randomly selected 
from the trajectory of MD simulation at 300K. Before alignment 
[Fig. 1(a)], superposition of all clusters translated to a coincident 
central point, with no rotation, shows that the neighboring atoms 
are distributed in uniform concentric spheres corresponding to the 
peaks of the Pair Correlation Function (PCF), with no resolved 
angular features. After the alignment, however, a clear 10-fold pat- 
tern emerges [Fig. 1(b)]. By applying a Gaussian smearing scheme to 
the 3D atomic distribution after the alignment, a three-dimensional 
(3D) spatially- resolved distribution function (SRDF) 23 D(r) can be 
defined in terms of the relevant 3D position vectors r rather than by 
simple distances as in a typical PCF analysis: 



-. mxn 



3/2 _«•(?-?,)* 



(1) 



where a is the Gaussian smearing parameter, r* is the position vector 
of atom i, m is the number of the clusters in the alignment and n is the 
number of atoms in each cluster. 

In Fig. 2, we present the iso-surfaces of the SRDFs for Cu centered 
clusters in Cu 64 . 5 Zr35 5 at five temperatures obtained from our align- 
ment scheme discussed above. The Cu-centered clusters show that 
ISRO is already well developed in the liquid at 1500K, which remark- 
ably, is —750 K above T g . As temperature decreases, the ordering 
becomes more pronounced as indicated by the increase in sharpness 
and intensity of the iso-surfaces. It is also clear that the icosahedral 
order gradually extends outward with decreasing temperature with 
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Figure 2 | Cu-centered SRDF at different temperatures (iso value : 0.081 A" 3 ). T m and T g are the melting point and the glass transition temperatures, 
respectively. Different colors are used here to indicate different coordination shells: red (first shell, icosahedron), blue (second shell, dodecahedraon), 
green (third shell, icosahedron), and white (fourth shell, rhombic triacontrahedron). The parameter a is chosen to be 4.0A" 2 to obtain optimal smoothness 
and resolution in the SRDF, but the resultant SRDF is not very sensitive to this parameter. The MRO at 300K can be seen to extend up to ~15A. 
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iso- surfaces clearly discernable out to as far as the fourth shell for the 
glass at 300 K. 

To quantify the SRO and MRO, we examine the intensities of 
individual peaks in the SRDF in each shell of the Cu-center cluster 
at four different temperatures. Fig. 3(b) shows the peak intensities at 
different temperatures averaged over all peaks in a given shell and 
normalized by D(r=0). The vertical bars in this figure indicate the 
variation of peak intensities within a given shell. Smaller variance at 
lower temperatures indicates stronger and more uniform icosahedral 
ordering. To obtain a characteristic decay length d for the icosahedral 
order at different temperatures, the directionally averaged peak 
intensities were fitted to an exponential decay function, l = e~ r l d , 
where r is the radius of a given shell in the SRDF (Fig. 3 and insert 
table). This provides a quantitative measure of the degree of ordering 
in this alloy. At lower temperatures, the decay length d increases due 
to increasing local order at longer lengths (Fig. 2). In the present 
analysis, the chemical identities of the atoms were not considered, 
but the extension of the analysis to obtain partial distributions is 
straightforward and is expected to yield insight into the chemical 
ordering in the system. 

Discussion 

Results from both scattering experiment and MD simulations both 
indicate a splitting in the second peak of the PCF in the Cu 6 4. 5 Zr 35 5 
system when cooling from the liquid to a glass 24 . This well-known 
splitting has been the subject of some speculation in the past 4 but 
there is as yet no firm evidence of the physical origin of this effect. 
The SRDF can be used to clearly identify the origin of this behavior 4 . 
Overlaying the atom positions for part of a Bergman triacontahedron 
cluster 25 along the 3 -fold direction over the SRDF, one can easily iden- 
tify the icosahedral-dodecahedral-icosahedral arrangement around 
the Cu center. Fig. 3(c) shows the Cu- centered partial PCF at different 
temperatures, where it is very clear that the second peak tends to split 
into two adjacent peaks (i.e. a right shoulder of the second peak 
develops) as temperature decreases. Our SRDF analysis reveals that 



the structural changes associated with the increasing second-peak 
shoulder corresponds to the development of strong ordering in the 
third coordination shell (green atoms in the Bergman model of 
Fig. 3(c)) below T g . This illustrates the utility of the cluster alignment 
scheme and the SRDF analysis, with which we can map the PCF into 
real space, thereby recovering the angular components of the radial 
distribution function. 

Finally, we also note that while the first shell 13 -atom icosahedral 
clusters are readily formed in the liquid at high temperatures, 
Bergman-type MRO only develops near T g . We examine the popu- 
lation and distribution of the Bergman-type clusters as the function 
of temperature. The atomic coordinates averaged over 200 successive 
MD snapshots at each temperature are used for the analysis in order 
to remove the effects from the thermal motions. The distribution of 
the Bergman-type clusters in the whole sample at 900, 700 and 300K 
are plotted in Fig. 4. At 900K, there are only some isolated Bergman- 
type clusters in the system (Fig. 4 (a)). Below Tg, a much higher 
density of bergman clusters is observed, forming a backbone network 
at 700 and 300K (Fig. 4 (b) and (c)). The clear emergence of a 
Bergman-type cluster network in the metallic glass suggests a pos- 
sible connection between the atomistic structures of metallic glasses 
and quasicrystals based on building blocks of Bergman clusters 25,26 . 
In the quasicrystal, the Bergman clusters are well oriented and have 
long range rotational order. While in metallic glasses, the orientation 
correlation between individual Bergman-type clusters is much 
weaker. The similarities and differences between the quasicrystals 
and metallic glasses may also give a clue to understanding why small 
variations in the chemical composition can alter the structure and 
formability of quasicrystal and metallic glass. For a given system, the 
chemical composition of the Bergman clusters should be similar and 
fixed more or less by the SRO in the system. However, additional 
small changes in the overall chemical composition can result in 
significant changes in the composition for the atoms in the glue 
region between the Bergman clusters and alter the spatial cluster- 
packing in favor or disfavor of phase ordering among the Bergman 
clusters. 




Figure 3 | A summary of the Cu-centered MRO at 300K. (a) SRDF (blue isosurface (0.081 A" 3 ), with local maxima shown as embedded spheres (Color 
represents different shells: red= 1st, blue=2nd, green = 3rd.); (b) radial decay of local maxima, averaged over each individual shell and normalized by the 
intensity at r=0, shown with exponential fit (The inset table shows the fitted decay lengths.); (c) correspondence between MRO structure and 
conventional Cu-Centered partial PCF: first shell (red spheres) to the first PCF peak, second shell (blue spheres) to the second PCF peak, and third shell 
(green spheres) to the developing third PCF peak. 
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Figure 4 | Development of Bergman triacontahedral clusters as the temperature is decreased, (a) 900K, 14 Bergman triacontahedral clusters which 
contain 12.1% of atoms; (b) 700K, 89 Bergman triacontahedral clusters which contain 46.9% of atoms; (c) 300K, 79 Bergman tiacontahedral clusters 
which contain 44.8% atoms in the whole sample. The Cu atoms are in blue, Zr in green and the center atoms of the Bergman triacontahedral clusters are 
in red. The population of the Bergman-type clusters is analyzed using the atomic coordinates averaged over 200 successive MD snapshots at each 
temperature to remove the effects from the thermal motions. The difference in the population of the Bergman-type clusters at 700 and 300K is partially 
due to the cooling from different initial configurations in the MD simulations. 



In summary, we provide a direct "order mining" tool that can be 
used to map out the SRDF, revealing and quantifying various aspects 
of MRO in noncrystalline systems which have remained poorly 
resolved with other analytic methods. We have applied the method 
to characterize the MRO that develops in Cu 6 4. 5 Zr 35 5 alloy as tem- 
perature decreases. We show that Bergman-type MRO and a back- 
bone network of interpenetrating Bergman-type clusters develop as 
the system approaches the glass transition. The discovery of Bergman 
triacontahedron MRO in the metallic glass sheds interesting light 
on the nature of the glass transition. In particular, it provides useful 
clues for understanding the kinetic pathway for slowing down the 
dynamics of the system and driving the system into non-ergodic 
behavior and the structural relation between metallic glasses and 
quasicrystals. Identification of these long range structural motifs also 
allow us to examine the influence of changing stoichiometry on the 
free energy competition between non- crystalline and various crys- 
talline motifs and hopefully understand the chemical influence on 
glass formation tendency. In addition to giving detailed description 
of the structural motifs developing in the system, the order para- 
meters obtained from such an "order mining" procedure can be 
incorporated into MD simulation (e.g., as collective variable in meta- 
dynamics) to overcome time-scale limitations in MD simulation of 
metallic liquid and glass systems. 

Methods 

In our alignment scheme, we characterize the local environment of each atom in the 
sample by an atomic cluster which includes the atom itself and its neighbors. The total 
number of such clusters is equal to the number of atoms in the sample. For binary 
alloys, we divide the clusters into subsets according to the chemical constituent of the 
center atom. The clusters in the same subset (according to the chemical constituent of 
the center atom) are aligned with each other to minimize the overall mean square 
distances between different clusters. The clusters are rigidly rotated and translated 
relatively to each other during the alignment. To carry out this alignment, a simulated 
annealing scheme implemented with molecular dynamics (MD) simulation is 
employed. The results of this alignment will reveal the common structural motif 
characterizing the local environment of a given chemical constituent. The use of the 
MD scheme permits the alignments of all the clusters in the subsets to be done 
simultaneously to obtain a consensus motif in an efficient manner. 

In the MD alignment simulations a modified Lenard-Jones type attractive potential 
is used to describe the attractive interaction between the atoms in different local 
clusters. 

W) = -°- 04 >< e >< (z^Ta) (2) 

\ r im,jn ' J 

where V im> j n is the attractive potential between atom m in cluster i and atom n in 
cluster; (i^j), and r inii j n is the distance between atom m in cluster i and atom n in 
cluster;'. The cutoff distance for the attractive interaction V im j n is set to be 1.0a, and 
the parameter in Equation (2) is set to 0.5a 2 . 



A strong harmonic potential is also introduced between atom pairs within 
the same cluster to keep the structure quasi-rigid during the collective -alignment 
process. 

U, n ,„=- l -Kx(r m ,„-rl n ) 2 (3) 

where U m>n is the strong harmonic potential between atom m and atom n, and the r m>n 
is the distance between atom m and atom n in the same cluster. r° m n is the initial 
distance between atom m and n before alignment. A force constant K=4.0X 10 4 s/a 2 
is large enough to guarantee that the geometric distortions of each cluster during the 
alignment are negligible. 

Although in principle all the clusters in a subset can be used in the collective- 
alignment, we find that 2000 local clusters randomly selected from each subset 
provide sufficient statistical sampling of the SRO in liquids and glasses. We also find 
that the alignments can be done more efficiently by first performing the alignment 
simulation on the small clusters involving only the first shell atoms, and then extend 
the alignment to further outer shells using the rotation matrices of the clusters 
determined from the first shell alignment. 
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